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Using 232 fb ^ of e^e collision data recorded by the BABAR experiment, we measure the rates 
of three-body Cabibbo-suppressed decays of the meson relative to the Cabibbo-favored decay: 

Iffr^-lt:"^) = (10.59 ± 0.06 ± 0.13) X 10-2 and ^gCy-^t;")' = (2-37 ± 0.03 ± 0.04) x 
where the errors are statistical and systematic, respectively. The precisions of these measurements 
are significantly better than those of the current world average values. 

PACS numbers: 13.25.Ft, 12.15.Hh, ll.30.Er 



Cabibbo-suppressed charm decays offer good labo- 
ratory for studying weak interactions as they provide 
a unique window on new physics affecting decay-rate 
dynamics and CP violation. The branching ratios of 
the singly Cabibbo-suppressed decays of D° meson are 
anomalous since the — > tt-tt^ branching fraction is 
observed to be suppressed relative to the _D" K~ 
by a factor of almost three, even though the phase space 
for the former is larger 0. The branching ratios of the 
three-body decays have larger uncertainties but do 
not appear to exhibit the same suppression. This mo- 
tivates the current study which measures the branching 
ratios of Z?" Tr^Tr+Tr^ and K^K'^n^ with respect to 
the Cabbibo-favored decay K^-k+t:^ Q. 

This analysis uses a data sample corresponding to an 
integrated luminosity of 232 fb^^ of e+e" collisions col- 
lected around ^/s w 10 GeV with the BABAR detector 
at the PEP-II asymmetric-energy storage rings. Track- 
ing of charged particles is provided by a five-layer sili- 
con vertex tracker (SVT) and a 40-layer drift chamber 
(DCH). Particle identification (PID) is provided by a 
likelihood-based algorithm which uses ionization energy 
loss in the DCH and SVT, and Cherenkov photons de- 
tected in a ring- imaging detector (DIRC). A large control 
sample of Z)*+ D^{K^ti^)tt^ events is used to eval- 
uate PID performance for kaons and pions from data. 
The average identification efficiencies and the misidenti- 
fication rates for pions (kaons) are 95% (90%) and 1% 
(3%) respectively. The typical separation between pi- 
ons and kaons varies from 8 standard deviations (cr) at 
momenta of 2 GcV/c to 2.5cr at 4 GcV/c. An electro- 
magnetic calorimeter (EMC) is used to identify electrons 
and photons. These systems are mounted inside a 1.5-T 
solenoidal magnet. The GEANT4 package is used to 
simulate the response of the detector with varying ac- 
celerator and detector conditions. Event reconstruction 
efficiency is obtained using Monte Carlo (MC) simulated 
events having production characteristics from the JET- 
SET fragmentation algorithm. Three-body decays 
are generated with uniform Dalitz plot (phase space) 
distributions. Electromagnetic radiation from final state 
charged particles (FSR) is modeled using the PHOTOS 
package g. 

To reduce combinatorial backgrounds, we reconstruct 
candidates in decays D*^ D^T^t i'^t is a soft, low 
momentum charged pion) with Z?" K^tt^tt'^, tt^tt^tt'^, 
and if-if+Tr", by selecting events with at least three 
charged tracks and a neutral pion. Photon candidates 



are reconstructed from calorimeter energy deposits above 
100 MeV, which are not matched to charged tracks. Neu- 
tral pions are reconstructed from pairs of photons with 
an invariant mass in the range 115-160 MeV/c^ and total 
energy in the laboratory system above 350 MeV. Charged 
kaon and pion and 7r+ candidate tracks are required to 
be within the fiducial volumes of the tracking and PID 
systems; they must have at least 20 hits in the DCH and 
transverse momenta greater than 0.1 GeV/c. Also, they 
must pass PID selection criteria. 

To form a Z)" candidate, two oppositely charged tracks 
and tt" are fit to a common vertex, constraining the 77 in- 
variant mass to the nominal tt" mass. The invariant mass 
of the candidate after the vertex fit is required to lie 
in the range 1.7-2.0 GeV/c^. To reduce high multiplicity 
events and combinatorial backgrounds, the momentum 
of the D° candidate in the event's center-of-mass frame 
(p*) is required to be greater than 2.77 GeV/c (this re- 
quirement also ensures that D° candidates from B decay 
are removed). The selected candidates after the above 
requirements are combined with the 7r+ track to form a 
D* + candidate. The D° and the 7r+ are constrained 
to originate from the coUison point; the resolution in 
Am, defined as the difference in invariant masses of the 
and candidates, is approximately 0.3 MeV/c^ 
for all three modes. Only those candidates are retained 
for which the vertex fit to the whole decay chain has a 
probability greater than 0.01 and Am is in the range 
144.9-146.1 MeV/c^ At this stage, approximately 3% of 
the events have multiple D*^ — !■ D^tt'^ candidates sat- 
isfying our selection criteria, due to misreconstruc- 
tion or a correctly reconstructed D° combining with a 
fake TTs- When there is more than one candidate in an 
event, we select only the candidate with the lowest ver- 
tex fit x^- Our selection procedures result in K~tt^tt^, 
TT-TT+TT^, and K~K^TT^ samples with purities 99%, 95% 
and 96%, respectively. 

The number of D'^ signal events in each decay mode 
is obtained by fitting the observed candidate mass 
distribution to the sum of signal and background com- 
ponents, where the latter has combinatorial contribu- 
tions and rcffcction contributions from real three-body 
D'^ decays where a kaon (pion) is misidentifed as a pion 
(kaon). The signal component is described by a sum of 
three Gaussians whose means and widths are allowed to 
vary. The combinatorial background is modeled by a lin- 
ear function. According to the MC simulation, a large 
fraction of the background consists of e~^e~ — > cc events. 
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with small contributions from processes e+e~ — * bb, uu, 
dd, ss. Reflected K~Tr'^Tr^ events peak in the lower 
(upper) sideband of rriT^-^+^o (m^^-^+^o). The levels 
of various background contributions in the Tr^Tr+Tr*^ and 
the K~ K'^tt'^ invariant mass distributions are shown in 
Fig. n The shapes of the K^tt^tt'^ reflections in the 
TT^TT^TT^ and K^K'^tt'^ invariant mass distributions are 
obtained from MC. The numbers of reflected events are 
found by making the A'^tt+tt" invariant mass distribu- 
tions for the 7r~7r^7r° and K~K~^'!t'^ samples and fitting 
them. Finally, maximum-likelihood fits are performed 
to extract the signal yields from the data for each of 
the three modes. The — > KqTt'^ decay is a Cabibbo- 



favored decay and a background for the D — > tt^tt+tt 
mode. The level of this contamination is obtained by 
peak in the m^-^+ distribution and the 



fitting the isTg 



number of Kgir^ events is subtracted from the 7r~7r+7r'' 
signal yield. The fitted candidate mass plots for the 
three modes are shown in Fig. [3 and the results of the 
fits are reported in Table 

The event reconstruction efficiency is obtained from 



Mode 


Number of 
signal events (S) 


Central value of 
D° mass (GeV/c^) 


Resolution 
(MeV/c^) 




505660 ± 750 


1.8646 ± 0.0002 


16.0 ± 0.5 




60426 ± 343 


1.8637 ± 0.0004 


17.4 ± 0.8 




10773 ± 122 


1.8649 ± 0.0004 


13.5 ± 1.0 



TABLE I: Number of observed signal events and the central 
value and resolution of the candidate mass distribution 
obtained from fit. The central value and the resolution are, 
respectively, the average mean and rms width of three Gaus- 
sians in the signal weighted by their fit-fractions. The errors 
are statistical only. 

MC. The reconstruction efficiency for each event is cal- 
culated as a function of its position in the Dalitz plot. 
That position is calculated using track momenta from a 
fit which constrains the h~h~^Tr'^ invariant mass to be the 
nominal D'^ mass, where 'h' is either a kaon or a pion. 
To correct for the differences in PID efficiency in data 
and MC, the ratio of these is determined for each track 
in bins of momentum and polar angle, and an event- 
by-event PID-correction factor is applied to each recon- 
structed event. The inverse of the calculated efficiency 
for each data point is taken as its weight. The average 
weight for each decay mode is computed by summing the 
weights of all events in the nominal signal regions (±3cr 
around the observed mean values of the mass distri- 
butions) and subtracting the efficiency-corrected event 
yields from sidebands (1.75-1.79 GeV/c^ and 1.95-1.99 
GeV/c^, spaced almost symmetrically around the nom- 
inal Z?" mass) to account for background events in the 
signal region. For the K~tt^tt'^ mode both sidebands are 
used for this purpose; for the 7r~7r+7r" (A'~A'+7r°) mode 
only the upper (lower) sideband is used because of the 
A'~7r+7r° reflection in the other sideband. The average 




1.85 1.9 1.95 
m iKTt^-nP) [GeV/c^] 



^ 10^ z 



10^ 




m (K"K"'7i;°) [GeV/c^ 

'° (top) and K 



K+TT° (bottom) 



FIG. 1: Simulated it' it ' n" 
invariant mass distributions. Signal events are shown as open 
histograms. Combinatorial and reflection backgrounds are 
shown by the light and dark shaded histograms, respectively. 
The signal region is delimited by the vertical lines. 



weights obtained from this method are verified to be un- 
biased. The average reconstruction weights for AT'tt+tt^, 
tt-tt+tt", and R-K+ir^ modes are 10.75 ± 0.02, 9.43 ± 
0.02, and 12.61 ± 0.05 respectively, where the uncertainty 
is due to MC statistics. 

The branching ratios are obtained from 



B{D° tt-tt+tt") _ N^-^+^o X W^-^+^o 



B{D^ A'-tt+ttO) Nk-^+^o X Wk~^+^o 



B{DO AT-tt+ttO) 



N 



K- 



X W, 



K- 



(1) 



(2) 



where N and W stand for the number of signal events 
detected and the average weight, respectively. 

The most important sources of systematic uncertain- 
ties in the branching ratios are reported in Table UTI The 
finite statistics of the MC samples used to obtain re- 
construction efficiencies contributes a small uncertainty. 
The uncertainty due to the Ato selection is estimated by 
repeating the analysis with different selection windows. 
The systematic uncertainty due to the background sub- 
traction procedure is obtained by repeating the analysis 
using cc Monte Carlo data and subtracting identifiable 
"true" background events in the signal region. The un- 
certainty caused by Dalitz plot binning effects in the ef- 
ficiency calculation is estimated by varying the bin-size. 
The effect of the modeling of the background probabil- 
ity distribution function on the signal yield is studied 
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1.75 



1.85 1.9 1.95 

m (K"K+7iO) [GeV/c^] 



FIG. 2: 



Fitted mass for the K tt n , n 



and 



K~K n data samples. Dots are data points and the sohd 
curves are the fit. The dot-dashed lines show the level of 
combinatorial background in each case. For the 7r~7r+7r° and 
the K~ K'^Tv° modes, the shaded region represents the total 
background. 



by repeating the fits to the candidate mass distri- 
butions with exponential and polynomial combinatorial 
background models. The systematic effect due to dif- 
ferences in the p* distribution in data and MC was de- 
termined by correcting the reconstruction efficiency ob- 
tained from MC by the the ratio of p* distributions in 
data and MC. Charged-particle identification studies in 
the data lead to small corrections applied to each track 
in the simulation. A large control sample of data and 
MC is studied separately to determine the residual PID 
uncertainties. Uncertainty due to potential differences 
in charged-particle tracking efficiencies in data and MC 
originating from an imprecise knowledge of different kaon 
and pion nuclear interaction cross sections and from the 
approximations used in our material model simulation, is 
conservatively assigned. 



As a consistency check, the analysis is performed sep- 



Si V p m fi i" 1 r* t; 




B(D"^K- K+Tv") 


S(D"-»K-7r+7r") 


S(D»^A'-7r + 7r») 


MC statistics 


0.27% 


0.47% 


Am selection 


0.30% 


0.90% 


Bg. Subtraction 


0.60% 


0.90% 


Efficiency binning 


0.11% 


0.24% 


Bg. PDF model 


0.16% 


0.13% 


p* difference 


0.24% 


0.02% 


PID 


0.77% 


0.84% 


Tracking 


0.60% 


0.60% 


Kg Removal 


0.07% 




Total 


1.25% 


1.73% 



TABLE II; Summary of systematic uncertainties. The total 
systematic uncertainty was obtained by adding the individual 
contributions in quadrature. 

arately for and D'^ events in different ranges of the 
D° candidate laboratory momenta to look for systematic 
variations as a function of charge or momentum outside 
the levels accounted for in the estimation of statistical 
and systematic uncertainties. The analysis is repeated for 
different data run periods and on the MC sample treated 
as data. As yet another cross-check, the branching ratios 
arc measured by directly fitting the efficiency-corrected 
histograms of the D*^ invariant mass distributions and 
then taking the ratio of the yields obtained from the fit. 
The results from all these cross-checks are consistent with 
the results of the main analysis. 

Using equations n and |21 we obtain the following re- 
sults for the branching ratios: 



B{D" tt-tt+tt") 



B{D" 



A'^tt+tt") 



= (10.59±0.06± 0.13) X 10^2^ (3) 



(2.37 ±0.03 ±0.04) x 10-^ (4) 



where the errors arc statistical and systematic, respec- 
tively. The previous most precise measurements for 
these branching ratios are (8.40 ± 3.11) x 10~^ and 
(0.95 ± 0.26) X IQ-^, respectively [13. We note that the 
second result differs significantly from the current world 
average value. As we consider events with any level of 
FSR as parts of our signals, the ratios we measure cor- 
respond to those of the so-called "bare" decay rates dis- 
cussed, for example, in Ref . • Using the world average 
value □ for the D° AT" 'TT'^TT^ branching fraction, we 
obtain, 

B{D° tt^tt+tt") = (1.493±0.008±0.018±0.053)xl0-^ 

(5) 

B{D" K-K+TT°) = (0.334±0.004±0.006±0.012)xl0-^ 

(6) 
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where the errors are statistical, systematic, and due to 
the uncertainty of B{D° K-Tr+n°). 

The decay rate for each process can be written as: 



r = J d^\M\^, 



(7) 



where F is the decay rate to a particular three-body fi- 
nal state, A4 is the decay matrix clement, and (f> is the 
phase space. Integrating over the Dalitz plot assuming a 
uniform phase space density, the above equation can be 
written as: 



{\M\' 



X $, 



(8) 



where (|A^p) is the average value of |A^p over the Dalitz 
plot and the three-body phase space, $ is proportional 
to the area of the Dalitz plot. For the three signal decays 
$ is in the ratio 7r"7r+7r° : K^n+n" : K^K+t:° = 5.05 
: 3.19 : 1.67. Combining the statistical and systematic 
errors, we find: 



{\M\^){DO ^ K-TT+TT°) 



(6.68 ±0.04 ±0.08) x 10"^ 

(9) 



= (4.53 ±0.06 ±0.08) x 10" 



(10) 



to date. The average squared matrix elements for both of 
the singly Cabibbo-suppressed decays are roughly a fac- 
tor of sin^ 9c smaller than that for the Cabibbo-favored 
decay and are therefore, in contrast to the corresponding 
two-body modes, consistent with the naive expectations. 
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<MW;^J™=,„8±0.14±0.21)xlO-.. 

(11) 

To the extent that the differences in the matrix ele- 
ments are only due to Cabibbo-suppression at the quark 
level, the ratios of the matrix elements squared for singly 
Cabibbo-suppressed decays to that for the Cabibbo- 
favored decay should be approximately sin^ 9c ~ 0.05 
and the ratio of the matrix elements squared for the two 
singly Cabibbo-suppressed decays should be unity. The 
deviations from this naive picture are less than 35% for 
these three-body decays. In contrast, the correspond- 
ing ratios may be calculated for the two-body decays 
K-TT+, and D° K' K+ . Us- 



ing the world average values for two-body branching ra- 
tios 0, the ratios of the matrix elements squared for 
two-body Cabibbo-suppressed decays, corresponding to 
Eas.lMTTl are, respectively, 0.034± 0.001, 0.111 ±0.002, 
and 3.53 ± 0.12. Thus the naive Cabibbo-suppression 
model works well for three-body decays but not so well 
for two-body decays. 

In summary, we have measured the ratios of the decay 
rates for the three-body singly Cabibbo-suppressed de- 
cays Tr^TT+TT*^ and K^K^n^ relative to that 
for the Cabibbo-favored decay K~Tr~^Tr^. This con- 
stitutes the most precise measurement for these channels 
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